Quantum model of classical mechanics: 
maximum entropy packets 

P. Hajicek 
Institute for Theoretical Physics 
University of Bern 
Sidlerstrasse 5, CH-3012 Bern, Switzerland 
hajicek@itp.unibe.ch 

April 2009 
PACS number: 03.65.Ta, 02.50.Tt 

Abstract 

In a previous paper, a statistical method of constructing quantum models 
of classical properties has been described. The present paper concludes the 
description by turning to classical mechanics. The quantum states that max- 
imize entropy for given averages and variances of coordinates and momenta 
are called ME packets. They generalize the Gaussian wave packets. A non- 
trivial extension of the partition-function method of probability calculus to 
quantum mechanics is given. Non-commutativity of quantum variables limits 
its usefulness. Still, the general form of the state operators of ME packets is 
obtained with its help. The diagonal representation of the operators is found. 
A general way of calculating averages that can replace the partition function 
method is described. Classical mechanics is reinterpreted as a statistical the- 
ory. Classical trajectories are replaced by classical ME packets. Quantum 
states approximate classical ones if the product of the coordinate and mo- 
mentum variances is much larger than Planck constant. Thus, ME packets 
with large variances follow their classical counterparts better than Gaussian 
wave packets. 



1 Introduction 



'The quantum origin of the classical' [TJ is a non trivial open problem of quantum 
theory: 'how to explain within quantum theory the classical appearance of our 
macroscopic world' [2]. The purpose of the present paper is to explain the classical 
properties as specific properties of quantum systems. The conceptual structure 
of quantum mechanics and with it the foundation of modern physics cannot be 
completely understood without such explanation. 

Every existing attempt in this direction starts with the assumption that the basic 
properties of individual quantum systems are single values of observables and that all 
other properties can be constructed or derived from these basic ones. Then, quantum 
mechanics does not admit any genuine realist interpretation. The definitive account 
is given by the Bub-Clifton-Goldstein theorem [3J. Only different kinds of apparent 
realism for different restricted sets of properties are possible. It is then difficult to 
explain how objective classical properties can emerge within quantum mechanics. 

In j3], we have initiated a very different approach. Our two main starting points 
are: 

1. Value of an observable o of an individual quantum system A is not a property 
of A alone, but of a composed system A + M, where M is an apparatus 
measuring o. The value of o measured by M is not determined (in general) 
before its measurement by M. It is created by the measurement process. We 
call, therefore, single values of observables measurable on A extrinsic properties 
of A. The whole existing quantum mechanics is practically only the theory of 
the extrinsic properties. 

2. As properties of a quantum system A, we allow also quantities that a) have 
values that may be more complex mathematical objects than just real numbers 
(such as sets, mappings between sets, etc.) and b) such a value need not be 
directly observable in a single measurement. 

Then, there are properties of quantum systems that can be viewed as determinate 
before their measurements without any other condition. They have been called in- 
trinsic, listed and classified into structural and conditional in [I]. Structural are 
those that are uniquely determined by the kind of quantum system (systems of the 
same kind are indistinguishable in the well-known strong and exclusively quantum- 
mechanical sense). Conditional are those that are uniquely determined by prepara- 
tions. In [4], a new realist interpretation of quantum mechanics has been described 
based on the intrinsic properties. 

Some of the intrinsic properties have been proposed as quantum models of clas- 
sical properties in [3]. In particular, two kinds of conditional properties have been 
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important: averages of quantum observables (including their variances) in a pre- 
pared state and the von-Neumann entropy of the state. Some explanatory remarks 
may be helpful. First, any preparation is defined by physical (objective) conditions. 
It need not be a process carried out by humans. Second, in the approach of [3], an 
average of an observable is not constructed or derived from the "more basic" single 
values of the observable. Averages are determined uniquely by preparation, single 
values are not. Not the single values determine the average but the average restrict 
possible single values. 

Third, entropy is often considered as a measure of observer ignorance and, there- 
fore, as a subjective concept. The "therefore" is fallacious. Entropy can generally 
be defined as follows. Let T be a complete set of mutually exclusive properties or 
states of a system A that is a measurable set and let the measure be ji{p) for p G T. 
Let physical condition C on A lead to restriction on possible properties or states 
of A such that the probability of the state p to occur under C is p(p). Then the 
entropy S(C) = -J r p(p)p(p). 

Such entropy can be a measure of ignorance in the following sense: everything 
we may know on A are some physical properties that therefore define some physical 
condition C . Then the above value of entropy gives the objective uncertainty on the 
properties or states of the system associated with condition C. Thus it is, in this 
case, if we know only C, simultaneously a measure of our ignorance on the system. 

Let us now briefly review the most popular approaches to the problem of classi- 
cality. At the present time, the problem does not seem to be solved in a satisfactory 
way, the shortcoming of the approaches being well known [5], [TU [15] . We mention 
them only fleetingly. First, the quantum decoherence theory [UI2] works only if cer- 
tain observables concerning both the environment and the quantum system cannot 
be measured (see the analysis in [51 [3] ) . The deep reason is that one works with val- 
ues of observables. Second, the theories based on coarse-grained operators [El [TJ [8] : 
the problem is the same as with the decoherence. For example, the Legget-Garg 
inequality [8] is a condition for the validity of the principle of macroscopic realism 
that works with values of observables. Third, the Coleman- Hepp theory [9j QUI [11] 
and its modifications [121 H3]: they are based on some particular theorems that 
hold for infinite systems but do not hold even approximately for finite ones (see the 
analysis in [TO]). 

The approach of [I] is free of these shortcomings. First, intrinsic properties are 
quantum properties of all quantum systems and there is no question about how 
they emerge in quantum mechanics. This avoids e.g. the artificial construction of 
classical propeties in the Coleman-Hepp approach. Second, they are considered as, 
and proved to be, objective in [3]. Hence, second, they could in principle serve as 
classical properties because they can satisfy the principle of classical realism. This 
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avoids the problems of both the quantum-decoherence and the coarse-grained theory 
that assume values of quantum observables to be basic properties. 

More specifically, [I] has conjectured that intrinsic averages and entropy for cer- 
tain macroscopic quantum systems can model all their classical properties. Clas- 
sical states of a macroscopic quantum system T have been defined as determined 
by averages Oi, - • • ,Ok of quantum operators ox, ■ ■ ■ , Ok that form a small subset 
of the algebra of observables of T. Finally, our modelling or construction of classi- 
cal properties is nothing but statistical physics. The statistical methods that were 
highlighted in [I] can only work if the following hypothesis (basic hypothesis of sta- 
tistical physics) is correct: The overwhelming part of macroscopic systems occur 
in quantum states that maximize entropy under the conditions of given averages 
Ox, - • - , Ok- This hypothesis is supported by observation and can be derived from 
quantum mechanics for a class of simplified models such as [TBI E3, EE] ■ 

The problem of quantum measurement (see, e.g., [SI EI3 El [IS] and references 
listed there) can be formulated within our approach as follows. For a measurement 
by a quantum apparatus M on a quantum system A, there must be an interaction 
between M and A as well as processes in M satisfying the conditions: (a) M changes 
its classical state as the result of the interaction, (b) the change depends on the initial 
quantum state of A, (c) average values defining different resulting classical states 
must differ by much more than the values of their variances, and (d) Born rule is 
fulfilled. The knowledge of what properties of M can be considered as classical is the 
first step. Still, to construct a model of interaction and processes in M satisfying 
conditions (a), (b), (c) and (d) remains a non trivial problem. Ref. [3] and the 
present paper do not offer a solution to the problem of measurement. 

A way of model construction for internal (thermodynamic) properties of macro- 
scopic quantum systems was described in detail in [I]. This, however, did not work 
for external (mechanical) properties of such systems. The aim of the present paper 
is to fill in this gap. 

The plan of the paper is as follows. In Sec. 2, classical mechanics is interpreted as a 
statistical theory. The existence of sharp trajectories is rejected so that all possible 
states of systems are described by fuzzy distribution functions. One choice for 
such distributions are the so-called maximum-entropy packets (ME-packet). These 
are states that maximize entropy for given averages and variances of coordinates 
and momenta. The method of partition function is used to calculate the general 
form of the distribution function. For a simple solvable example, the dynamical 
equations for the averages and variances are obtained. The example shows how the 
equations of motion are reinterpreted in our theory. For general potentials, we use 
an approximative method: step-by-step calculation of the higher and higher time 
derivatives of coordinates and momenta. This will later be compared with quantum 
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ME packets. 

Turning to quantum mechanics, we apply the maximum-entropy principle in an 
analogous way in Sec. 3. The averages and variances are taken over from the classi- 
cal states that are to be modelled. A straightforward generalization of the partition 
function method is now complicated by the non-commutativity of the coordinates 
and momenta. We can show that only the first derivatives of the logarithm of par- 
tition function have the usual meaning. This is, however, sufficient for calculating 
the state operators for all ME packets. We find the diagonal representation of the 
state operator in Sec. 3.2 and obtain with its help the general form of the partition 
function and the state operator itself. It turns out that Gaussian wave packets are 
special case of ME packets, namely those with zero entropy and minimal uncer- 
tainty. The diagonal representation gives us also a powerful method to calculate 
averages of higher moments. In fact, what has been done in Sees. 3.1 and 3.2 is a 
non-trivial extension of the partition-function method of the probability calculus as 
described, e.g., in [20] to quantum mechanics that might also be of some interest for 
mathematicians. 

In Sec. 3.3, the equations of motion are calculated in analogy to the classical 
case. We find that the quantum corrections to the classical equations come only 
from high powers of q in the expansion of the potential or in high powers of t in 
the expansion of the time-dependent averages. Also, these corrections are of the 
second order in H. These results show that our quantum models follows classical 
trajectories very closely. The nature of classical limit is studied in Sec. 4. The result, 
which may seem surprising, is that it is the limit of large variances, not small. Thus, 
quantum ME packets with large variances follow their classical counterparts better 
than Gaussian wave packets. Of course, the way we measure the size of the variances 
is important here. The variances that are large with respect to this measure can 
still be sufficiently small to agree with observations. Finally, Sec. 5 concludes the 
paper by summarizing the main ideas and the most important results. 

2 Statistical form of classical mechanics 

Let us start with the warning that the topic of this section has nothing to do with 
what is usually called 'statistical mechanics'. 

If one is going to model classical mechanics then what are the properties that one 
would like to reproduce? The most conspicuous property from the point of view of 
quantum mechanics appears to be the sharpness of mechanical trajectories in the 
phase space because quantum mechanics denies the existence of such trajectories. 
This leads most researchers to aim at quantum states the phase-space picture of 
which is as sharp as possible. That are states with minimum uncertainty allowed 
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by quantum mechanics. For one degree of freedom, described by coordinate q and 
momentum p, the uncertainty is given by the quantity 



It is well known that minimum uncertainty allowed by quantum mechanics is v — 1. 

The states with v = 1 are, however, very special states. First, they must be pure 
states such as Gaussian wave packets or coherent states. Such states are very difficult 
to prepare unlike the usual states of macroscopic systems described by classical 
mechanics. They are also prone to strong distortion by measurements. Moreover, 
as pure states, they can be linearly superposed. This is another peculiarity that is 
never observed for states of systems of classical mechanics. Hence, trying to get a 
trajectory as sharp as possible leads to the loss of other desirable properties. 

Moreover, observations within classical mechanics admit the notion that the 
sharpness of phase-space trajectories is only a mathematical and methodical fea- 
ture of classical mechanics. It may be just an idealisation, a limit in which things 
become mathematically simpler. We can use it in calculations which, however, must 
also take into account the necessary non-zero variances of real observations. In- 
deed, such observations are generally afflicted with uncertainties v ^> 1. Hence, if 
we want to compare the predictions of our quantum models with observations of 
classical mechanics, we are forced to compare states that are fuzzy in both theories. 

One idea of the present paper is to consider states with given averages and vari- 
ances of the coordinates and momenta and leave everything else as fuzzy as possible. 
To calculate the corresponding probability distributions in classical, and the state 
operators in quantum mechanics, we shall, therefore, apply the maximum entropy 
principle. This is a general principle in mathematical theory of probabilities (see 
|20j) and it should not be confused with the well-known thermodynamic law. The 
resulting states are called maximum- entropy packets, ME-packets. The averages of 
coordinates and momenta take over the role of coordinate and momenta in classical 
mechanics. In any case the averages represent measurable aspects of these variables. 
The dynamical evolution of variances is an important indicator of the applicability 
of the model one is working with. It determines the time intervals within which 
reasonable predictions are possible. 

Consider a three-body system that is to model the Sun, Earth and Jupiter. It 
turns out that generic trajectories starting as near to each other as, say, the di- 
mension of the irregularities of the Earth surface will diverge from each other by 
dimensions of the Earth-Sun distance after the time of only about ten million years. 



2AqAp 




v = 



h 



where Aa is the variance of quantity a, 



Aa = vV) - (a) 2 . 



(2) 
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This seems to contradict the four billion years of relatively stable Earth motion 
around the Sun that is born out by observations. The only way out is the existence 
of a few special trajectories that are much stabler than the generic ones and the 
fact that bodies following an unstable trajectory have long ago fallen into the Sun 
or have been ejected from the solar system. 

An important question is that on the ontological status of ME-packets and on 
the nature of the limit in which trajectories become sharp. The usual standpoint is 
that any mechanical system always objectively is in a state of a completely sharp 
trajectory. Any more fuzzy state is only the result of our incomplete knowledge. 
Thus, the fuzzy states are not considered themselves as real. Here, we take the 
opposite standpoint. For us, a state to be real, it must be determined by objective 
initial conditions. A simple example is a gun in a position that is fixed in a re- 
producible way and that shoot bullets using cartridges of a given provenance. The 
state of each individual shot is defined by the conditions and is the same for all shots 
even if observations may have different results for different shots. A finer analysis is 
possible only as long as new initial conditions are specified that determine a subset 
of individual shots. In the theoretical description of a state, we can make the limit 
of AQ — > 0, AP — > 0. This is considered as a non-existing, but practically useful 
idealization. 

To limit ourselves just to given averages and variances of coordinates and mo- 
menta is a great simplification that enables us to obtain interesting results easily. 
Some further discussion on quantum modelling of classical properties is in order. On 
the one hand, for internal degrees of freedom, the usual thermodynamic methods 
give small relative variances as a consequence of the state coordinates being exten- 
sive, the entropy being maximal and the system being macroscopic. This does not 
work for external (mechanical) properties. The difference is due to the simple fact 
that the internal degrees of freedom are not accessible to manipulation and have 
small variances spontaneously. The external degrees of freedom are accessible to 
manipulations and it is easy to prepare states with small as well as large variances. 
There is no objective need that the variances are small spontaneously. The idea that 
really existing mechanical states must always have only small variances is caused 
by a purely theoretical notion that all real mechanical systems have an absolutely 
sharp phase space trajectory and this notion is clearly false. The only problem is 
that it has become a part of our subconscious psychology. 

On the other hand, if quantum ME packets are to be quantum models of the 
classical ME packets with the same averages and variances then this is a more 
general situation than that considered in Ref . [I] . Quantum ME packet is a classical 
state in the sense of Ref. [I] if it has small variances. Only then, the average values 
are directly observable on individual systems. 
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2.1 Classical ME-packets 

Let us first consider systems with one degree of freedom. The generalization to 
any number is easy. Let the coordinate be q and the momentum p. A state is a 
distribution function p(q, p) on the phase space spanned by q and p. The function 
p(q,p) is dimension-free and normalized by 

dq dp 

P = 1 , 

v 

where v is an auxiliary phase-space volume to make the integration dimension-free. 
The entropy of p(q,p) can be defined by 

f dqdp 
b := — pmp . 



The value of entropy will depend on v but the most of other results will not. Classical 
mechanics does not offer any idea of how to fix v. We shall get its value from quantum 
mechanics. 

Let us define: ME-packet is the distribution function p that maximizes the en- 
tropy subjected to the conditions: 

(q) = Q, (q 2 ) = AQ 2 + Q 2 , (3) 

and 

(p) = P, (p 2 ) = AP 2 + P 2 , (4) 

where Q, P, AQ and AP are given values of averages and variances of q and p. We 
have used the abbreviation 

dqdp 



(x) = / x . 

J v 

The explicit form of p can be found using the partition-function method as de- 
scribed e.g. in [20]. The variational principle yields 



p = ~w\ — \ — \ — rr ex p(-^i<? - - hq 2 - Kp 2 ) , (5) 

^ [M, -A2, A3, A4J 

where 

^l^R. eX p(-Aig - \ 2 p - \ 3 q 2 - A 4 p 2 ) , 

v 

and Ai, A2, A3 and A4 are the Lagrange multipliers. Hence, the partition function 
for classical ME-packets is given by 
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The expressions for Ai, A2, A3 and A4 in terms of Q, P, AQ and AP can be obtained 
by solving the equations 



and 



The result is: 



and 



dAi ^ ' <9A 3 



dlnZ dlnZ A 9 9 

aA 2 aA 4 



A - P A ' 



AP 2 ' * 2AP 2 

Substituting this into Eq. ([2]), we obtain the distribution function of a one-dimensional 
ME packet. The generalization to any number of dimensions is trivial. 

Theorem 1 The distribution function of the ME-packet for a system with given 
averages and variances Qi, ■ ■ ■ , Q n , AQi, ■ • • , AQ n of coordinates and Pi, ■ • • , P n , 
APi, • ■ • , AP n of momenta, is 



n f j 



(?*; - Qk) 2 (pk - Pk) 



2AQ 2 k 2AP 2 



(9) 



We observe that all averages obtained from p are independent of v and that the result 
is a Gaussian distribution in agreement with Jaynes' conjecture that the maximum 
entropy principle gives the Gaussian distribution if the only conditions are fixed 
values of the first two moments. 

As AQ and AP approach zero, p becomes a delta- function and the state becomes 
sharp. For some quantities, this limit is sensible for others it is not. In particular, 
the entropy, which can easily be calculated, 

g=1+ln 2^AgAP 

V 

diverges to —00. This is due to a general difficulty in giving a definition of entropy 
for a continuous system that would be satisfactory in every respect. What one could 
do is to divide the phase space into cells of volume v so that AQAP could not be 
chosen smaller than v. Then, the limit AQAP — > v of entropy would make more 
sense. 

The average of any monomial of the form q k p l q 2m p 2n can De calculated with the 
help of partition-function method as follows: 

(a k n l a 2m v 2n ) = ^ — (10) 
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where N = k + I + 2m + 2n, Z is given by Eq. and the values ((Tj) and (jHJ) must 
be substituted for the Lagrange multipliers after the derivatives are taken. 

Observe that this enables to calculate the average of a monomial in several dif- 
ferent ways. Each of these ways, however, leads to the same result due the the 
identities 

d 2 Z dZ d 2 Z dZ 
which are satisfied by the partition function. 
2.2 Equations of motion 

Let us assume that the Hamiltonian of our system has the form 

* = £ + (id 

where m is the mass and V(q) the potential function. The equations of motion are 

<1 = {(1,H}, p={p,H}. 
Inserting ( TTTI) for H, we obtain 

. P . dV 

q = — , p = -~r ■ (12) 

m aq 

The general solution to these equations can be written in the form 

q(t) = q(t;q,p) , p(t) = p(t;q,p) , (13) 

where 

q(0; q,p) = q, p(0; q,p)=p, (14) 

q and p being arbitrary initial values. We obtain the equations of motion for the 
averages and variances: 

Q(t) = (q(t;q,p)) , AQ(t) = a/ ((q(t; q,p) - Q(t)) 2 ) (15) 

and 

P(t) = (p(t;q,p)) , AP(t) = ^((p(t;q lP )-P(t)) 2 } . (16) 

In general, Q(t) and P(t) will depend not only on Q and P, but also on AQ and 
AP. 

Let us consider the special case of at most quadratic potential: 

V(q) = V + V iq +^V 2 q 2 , (17) 
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where Vk are constants with suitable dimensions. If V\ = V 2 = 0, we have a free 
particle, if V 2 = 0, it is a particle in a homogeneous force field and if V 2 7^ 0, it is an 
harmonic or anti-harmonic oscillator. 

In this case, the general solution has the form 

q(t) = fo(t)+qfi(t)+pf 2 (t) , (18) 
p(t) = g {t)+qgi{t)+pg 2 (t) , (19) 

where / (0) = / 2 (0) = g Q (0) = ^(0) = and /x(0) = g 2 (0) = 1. If V 2 ^ 0, the 
functions are 

Vi 1 
fo(t) = - — (1 - cosut) , fi(t) = coswt , f 2 (t) = - smut , (20) 
v 2 4 

W 

9o(t) = -£—smut , g x (t) = -£smut , g 2 (t) = cos cut , (21) 
V 2 

where 



2 , u 




m 



Only for V 2 > 0, the functions remain bounded. If V 2 = 0, we obtain 

/o(*) = -^. fi(t) = l, /„(*) = -, (22) 
2m m 

go (t) = -V x t , ^(i) = , g 2 {t) = 1 . (23) 

The equations for averages and variances resulting from Eqs. (TT3"|) . (EJ) and (J4j) 

are 

Q(t) = fo(t)+Qfx(t)+Pf 2 (t) , (24) 

and 

AQ 2 (t) + Q 2 (i) = f {t) + (AQ 2 + Q 2 )/ 2 (t) + (AP 2 + P 2 )/ 2 (t) 

+ 2Qf (t)f 1 (t) + 2Pf (t)f 2 (t) + 2{qp)f x {t)f 2 {t) . (25) 

For the last term, we have from Eq. (11 01) 

1 d 2 Z 



(qp) 



Z d\id\ 2 

Using Eqs. ©, © and ©, we obtain from Eq. p5|) 



aq(o = jfm&Q 2 + fimp 2 ■ (26) 



10 



Similarly, 

P(t) = g Q (t) + Qgi(t) + Pg 2 (t) , (27) 
AP(t) = yj f]{t)AQ 2 + g%(t)AP 2 . (28) 

We observe: if functions fi(t), f2(t), gi(t) and g^it) remain bounded, the variances 
also remain bounded and the predictions are possible in arbitrary long intervals of 
time. Otherwise, there will always be only limited time intervals in which the theory 
can make predictions. 

In the case of general potential, the functions (TT3T) can be expanded in products 
of powers of q and p, and the averages of these products will contain powers of the 
variances. However, as one easily sees form formula (fTUl) and (jSJ), 

(gV) = Q k p l + XAQ + YAP , 

where X and Y are bounded functions. It follows that the dynamical equations 
for averages coincide, in the limit AQ — > 0, AP — > 0, with the exact dynamical 
equations for q and p. It is an idealisation that we consider as not realistic, even in 
principle, but that may still be useful for calculations. 

Let us expand a general potential function in powers of q, 

v ® = E k\ Vkqk > ( 29 ) 

k=0 

where are constants of appropriate dimensions. The Hamilton equations can be 
used to calculate all time derivatives at t — 0. First, we have 

-dt =[q ' H} = m- 

This equation can be used to calculate all derivatives of q in terms of those of p: 

dt n mdt n - 1 ' ^ ' 

A simple iterative procedure gives us further time derivatives of p: 

d P T/ T/ V 3 2 V i 3 , /oi^ 

"77 = ~Vi - V 2 q ~—q- —q a + r 5 , 31 
at 2 b 

d 2 P V 2 V 3 V, 2 

"TTT = P qP - 7^—1 P + r 5 , (32) 

dt z m m 2m 

d 3 P v 3 2 7 4 2j to , kv 3 + v^ 2 , 3y 2 y 3 + 2 

-77T = 2^P + + 9 + o 9 

at J m^ m m 2m 

+ 4V^+3V£ 3 + 5VM 4 + + ^ 

6m 12m 12m ' 
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and 



d*p V 4 , + 3VxV 4 + W 2 V 3 5V? + 8V 2 V A , 

171 = 3P + 5 P + 2 qp + ^ g P 

at 4 m * m z m A 2m z 

v 3 v 4 3 , wi 4 



+ 3^g^+^ 4 p + r 5 , (34) 



m 



where is the rest term that is due to all powers in (129|) that are not smaller than 
k (the rests symbolize different expressions in different equations). The purpose of 
having all time derivatives up to the fourth order is to show later that it is the highest 
order in which no quantum corrections appear in the equations for the averages. 

Taking the average of both sides of Eqs. (|5I)) - (|33)I . and using Eq. (TJO]), 
we obtain 

£ = - V, - V 2 Q - - - Yl±m^ + ,, , (35) 

'£L = -Yip + Yi Q p + Yl q ^p + Y± PAQ i + rb , (36 ) 

dt z m m 2m 2m 

d " P _ ^3„ 2 Vj„ n2 , W 2 , VW + lf n , 3^3 + W 4 n2 

ctr d rrr m/ m m 2m 

6m 12m 12m \m/ / 

+ / 3W 3 + W 4 + W 2 V 4 + Wj Q + 5^ g2 + 5^ aq2 
\ 2m 2m 2m 4m 

c;t/ 2 W 2 \ 
+ TT^Q 3 + X^QAQ 2 AQ 2 + r 5 , (37) 
bm 4m / 

and 



= _ii P 3 + WV+lf p + 3V 1 V 4 + 5V 2 V 3 
dt 4 m 3 m 2 m 2 



+ 5V?+^ Q2p + ^ + 3^ _ 3^ 
2m" 1 m 2 4m z m * 

5K 2 + 8V 2 \/ 4n 9V 3 7 4/ , n 9F 4 2 „ 9V 4 2 



+ %' 2 * P + ~^QP + ^~ 2 Q 2 P + ~r- 2 P^Q 2 + r 5 • (38) 

\ 2m z m z 2m z Am z / 

We can see, that the limit AQ — > 0, AP — > in Eqs. (I3"51) - (13"8"l) lead to equations 
that coincide with Eqs. ([3T]) - (jM|l if Q — ► q, P — > p as promised. 



3 Quantum ME-packets 

Let us now turn to quantum mechanics and try to solve an analogous problem. Let 
a system with one degree of freedom be described by the operators q and p and let 



12 



us look for a state p, a normalized, 



Tip = 1 , 



self-adjoint positive operator, that maximizes von Neumann entropy 



S 



Tr(plnp) 



(39) 



under the conditions 



Tr(pg) = Q , 



Tr(pg 2 ) = Q 2 + AQ 2 



(40) 



Tr (pp) = P 



Tr(pp 2 ) = P 2 + AP 2 



(41) 



where Q, P, AQ and AP are given numbers. The states that satisfy these conditions 
are called quantum ME-packets. 

3.1 Calculation of the state operator 

To solve the mathematical problem, we use the method of Lagrange multipliers as 
in the classical case. Thus, the following equation results: 

dS - A dTrp - AidTr(pg) - A 2 dTr(pp) - \ 3 dTr(pq 2 ) - \ 4 dTr{pp 2 ) = . (42) 

The differentials of the terms that are linear in p are simple to calculate: 



Although not all elements of the matrix dp mn are independent (it is a hermitian ma- 
trix), we can proceed as if they were because the matrix x nm is to be also hermitian. 
The only problem is to calculate dS. We have the following 

Lemma 1 



Proof Let M be a unitary matrix that diagonalizes p, 

M pM = R , 

where R is a diagonal matrix with elements R n . Then S — — ^2 n R n \nR n . Correc- 
tion to R n if p i— > p + dp can be calculated by the first-order formula of the stationary 
perturbation theory. This theory is usually applied to Hamiltonians but it holds for 




mn 




(43) 



mn 
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any perturbed hermitian operator. Moreover, the formula is exact for infinitesimal 
perturbations. Thus, 

Rn >-> Rn + M L M lndpkl ■ 

kl 

In this way, we obtain 



dS = - ( Rn + M L M lndpkl J 

n \ kl / 

RnU + ^J2 M rn M sndp r )j 

E 



/ J Rn m Rn 



In R n M\ n M ln dp kl + Mt n M ln dp, 

kl kl 

^ [hi + (hip)fc*]] dp 



kl 

Q.E.D. 

With the help of Lemma 1, Eq. (|4"2l) becomes 

Tr [(1 + In p - A - Xiq - A 2 p - A 3 g 2 - A 4 p 2 )(ip] = 
so that we have 

p = exp(-A - 1 - Aig - \ 2 p - A 3 g 2 - A 4 p 2 ) . (44) 
The first two terms in the exponent determine the normalization constant 

e -Ao-l 

because they commute with the rest of the exponent and are independent of the 
dynamical variables. Taking the trace of Eq. (j44j) . we obtain 

1 



-Ao-l 



Z (Ai, A2, A3, A4) 

where Z is the partition function, 

Z(X X , A 2 , A 3 , A 4 ) = Tr[exp(-Aig - \ 2 p - A 3 g 2 - A 4 p 2 )] . 
Thus, the state operator has the form 

^ [M, ^2, A3, A4J 

At this stage, the quantum theory begins to differ from the classical one. It turns 
out that, for the case of non-commuting operators in the exponent of the partition 
function, formula (jTOj) is not valid in general. We can only show that it holds for 
the first derivatives. To this aim, we prove the following 



(45) 



(46) 
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Lemma 2 Let A and B be hermitean matrices. Then 

Tr[exp(A + BX)} = Tr[B exp(A + BX)) . (47) 

dX 

Proof We express the exponential function as a series and then use the invariance 
of trace with respect to any cyclic permutation of its argument. 



dTr[exp(A + BX)] = ^ ^ Tr K^ + BX) n ] 

n 

^{A + BX) k - 1 B{A + BX) 



n=0 



n=0 



n—k 

fc=l 



dX 

]T — } Yl Tr \- B ( A + BX ) n ^] dX = Tl l B e M A + BX)]dX 



n. 

n=0 k=l 

Q.E.D. The proof of Lemma 2 shows why formula (TTUT) is not valid for higher deriva- 
tives than the first in the quantum case: the operator B does not commute with 
A + BX and cannot be shifted from its position to the first position in product 

(A + BX) k B(A + BX) 1 . 

For the first derivative, it can be brought there by a suitable cyclic permutation. 
However, each commutator [B, (A + BX)] is proportional to h. Hence, formula ffTUj) 
with higher derivatives is the leading term in the expansion of averages in powers of 

h. 

Together with Eq. (145]) . Lemma 2 implies the formulae: 

^ = -0, ^=-«*-AQ* (48) 

oX\ oX 3 

and 

dlnZ dhiZ 9.9 , s 

"aA7 = ' Hxr = (49) 

The values of the multipliers can be calculated from Eqs. (j4~8l) and (jUJ), if the form 
of the partition function is known. 

Variational methods can find locally extremal values that are not necessarily 
maxima. We can however prove that our state operator maximizes entropy. The 
proof is based on the generalized Gibbs' inequality, 

Tr (p In p — p In a) > 

for all pairs {p,u} of state operators (for proof of the inequality, see [6], P. 264). 
The proof of maximality is then analogous to the 'classical' proof (see, e.g., [20], P. 
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357). The first proof of maximality in the quantum case was given by von Neumann 

I2H- 

The state operator (146]) can be inserted in the formula (139]) to give the value of 
the maximal entropy, 

S = InZ + Ai(g) + \ 2 ( P ) + A 3 (g 2 > + A 4 (p 2 ) . (50) 

This, together with Eqs. (HHl) and (l4"9"|) can be considered as the Legendre transfor- 
mation from the function lnZ(Ai, A 2 , A 3 , A 4 ) to the function S((q), (p), (q 2 ), (p 2 ))- 



3.2 Diagonal representation 

The exponent in Eq. (1461) can be written in the form 

A + 3_- 2 ^x lK , (51, 

where 

*-i^(« + ^) , + 5^(' + ^)'- « K » 

This is an operator acting on the Hilbert space of our system. K has the form of 
the HamiltoniaE0 of a harmonic oscillator with the coordinate U and momentum W 

U = q + ^k' W = P+ ^t 4 ' (53) 

that satisfy the commutation relation [U, W] = ih. The oscillator has mass M and 
frequency Q, 

M=W^, Vl = l. (54) 
V A 4 

The normalized eigenstates \k) of the operator form a basis in the Hilbert space 
of our system defining the so-called diagonal representation and its eigenvalues are 
h/2 + hk. As usual, we introduce operator A such that 



u = 




(55) 


w = 




(56) 


K = 


£(AlA + AAi)), 


(57) 


A\k) = 


Vk\k-i) , 


(58) 


A*\k) = 


Vk+ l\k + 1) . 


(59) 



^^The operator K must not be confused with the Hamiltonian H of our system, which can be 
arbitrary. 
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To calculate Z in the diagonal representation is easy: 

00 / \2 \2 \ 



/ ^2 ^2 \ 00 

eXP V 4A~ + 4A" ^ ) J2 e M-2ft^f>^k) . 

V 3 4 J k=(j 



Hence, the partition function for the quantum ME-packets is 

ex P (Ss + Si 



Z = \ - J- . (60) 

2 sinh(^v^3^4) 

Now, we can express the Lagrange multipliers in terms of the averages and vari- 
ances. Eqs. (HHD and §M) yield 

Al = -A^2 ln U3I' A2 = -Ai^2 ln ^I' (61) 
and , 

A3= W2 ln ^I' A4= 2Ai^2 ln ^I' (62) 
where v is defined by Eq. (CQ) 

From Eq. ( !50l) . ( !6~TI) and f |62l) . we obtain the entropy: 

5 = _l n 2 + ^±lln(z/+l)-^yiln(z/-l) . (63) 

Thus, 5* depends on Q, P, AQ, AP only via v. We have 

dS 1 z/ + l n 

— = - In > , 

2 z/ - 1 

so that 5 1 is an increasing function of v. Near v — 1, 

lm>-l) . 

Asymptotically (^ — > 00), 

S« lni/ + l -ln2 . 

In the classical region, v 3> 1, S « InzA 

It is clear that the choice of Q and P cannot influence the entropy. The in- 
dependence of S from Q and P does not contradict the Legendre transformation 
properties. Indeed, usually, one would have 

dS _ 

ag" Al ' 
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but here 



H = X 1 + 2A 3 g 
oQ 



which is zero. 

The state operator can also be expressed in terms of the averages and variances. 
The trivial generalization to n degrees of freedom is 

Theorem 2 The state operator of the ME-packet of a system with given averages 
and variances Qi, ■ ■ ■ , Q n , AQi, ■ ■ ■ , AQ n of coordinates and Pi, ■ ■ ■ , P n , AP 1; ■ • • , 
AP n of momenta, is 



P 



n 

k=l 



l4 



exp 



1 , Vk 
-— In ■ 



ft v k -\ 



where 



and 



K 1APk f 
Ku — (Qk 

2 AQ k Kqk 



-Qk) 2 + \^(Pk-Pk) 2 

2AP k AQ k 
ft 



(64) 



(65) 



(66) 



Strictly speaking, the state operator ( 164"1) is not a Gaussian distribution. Thus, it 
seems to be either a counterexample to, or a generalization of, Jaynes statement 
that the Gaussian distribution is the only distribution that maximizes entropy for 
given values of the first two moments [20J. 
In the diagonal representation, we have 



p = J2R k \k)(k\ 



(67) 



k=0 



We easily obtain for R k that 



Pi 



jy- l) k 
{u + l) fc+r 



(68) 



Hence, 



lim R k = 5 k0 



and the state p becomes |0)(0|. In general, states \k) depend on v. The state vector 
|0) in the g-representation expressed as a function of Q, P, AQ and v is given by 



i){q) 



1 



tt2AQ 2 



1/4 



exp 



v 



AAQ 



iPq 
~ft~ 



(69) 



This is a Gaussian wave packet that corresponds to different values of variances 
than the original ME packet but has the minimal uncertainty. For v — > 1, it remains 
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regular and the projector |0) (0| becomes the state operator of the original ME packet. 
Hence, Gaussian wave packets are special cases of ME-packets. 

The diagonal representation offers a method for calculating averages of coordi- 
nates and momenta products that replaces the partition function way. Let us denote 
such a product X. We have 

oo 

(X) = ^2R k (k\X\k) . (70) 

k=0 

To calculate (k\X\k), we use Eqs. (155]), (I3B1), ([55]) . (1511). AST]) and (fS2]l to obtain 
q = Q + ^l(A + A^ , p = P-i^L(A-Ai). 



By substituting these relations to X and using the commutation relations [A, A*] = 
1, we obtain 

X = P(N) + Q(A, A*) , 

where N = A^A and where, in each monomial of the polynomial Q, the number of 
A-factors is different from the number of A^-factors. Thus, 

(k\X\k) = V(k) . 

In Eq. (170|) . there are, therefore, sums 



J2 knR k 



k=0 



With Eq. ( 1681) . this becomes 



E 2 
k Rk = ; T-^n ; 
V + 1 

k=0 



where 



We easily obtain 



oo / _ 

k=0 v 



,u 2 -ld\ n u + l 

J- II 



2 di/ 

The desired average value is then given by 

The calculation of the polynomial P for a given X and the evaluation of the right- 
hand side of Eq. ([71]) are the two steps of the promised method. 
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3.3 Equations of motion 

Let the Hamiltonian of our system be H and the unitary evolution group U(t). The 
dynamics in the Schrodinger picture leads to the time dependence of p: 

p{t) = U{t)pU\t) . 

Substituting for p from Eq. ( jMj) and using a well-known property of exponential 
function, we obtain 



p(t) = -^exp(~ln^U(t)KU\t^ . (72) 



v 2 



In the Heisenberg picture, p remains constant, while q and p are time dependent 
and satisfy the equations 

ih^ = [q,H], ihf t = [p,H}. (73) 

They are solved by 

q(t) = U\t)qU(t) , p(t) = U\t)pU{t) , 

where q and p are the initial operators, q = q(0) and p = p(0). The resulting 
operators can be written in the form of operator functions analogous to classical 
expressions ([TBI so that Eqs. (}T5]) and f|T6|) can again be used. 

The example with potential function (fT7j) is solvable in quantum theory, too, and 
we can use it for comparison with the classical dynamics as well as for a better 
understanding of the ME-packet dynamics. Eqs. (173|) have then the solutions given 
by (HBD and JTSJ with functions f n {t) and g n {t) given by fl2DD and (EJJ) or ([22]) and 
(1231) . The calculation of the averages and variances is analogous to the classical one 
and we obtain Eqs. (I24p and (I25p again with the difference that the term 2(qp) on 
the right hand side of (12"oT) is now replaced by (qp + pq) . 

To calculate (qp + pq), we use the method introduced in the previous section. We 
have 

qp + pq = 2QP + 2^3-(A + tf) - 2z^£(A - A*) - 2i^^(A 2 - A^ 2 ) . 



hence, V = 2QP, and 

(qp + pq) = 2QP . 

The result is again Eq. ( |26i) . Similarly for p, the results are given by Eqs. ( 1271) and 



We have shown that the averages and variances of quantum ME-packets have 
exactly the same time evolution as those of classical ME-packets in the special case 
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of at-most-quadratic potentials. From formulae (126]) and (I28I) we can also see an 
interesting fact. On the one hand, both variances must increase near t = 0. On the 
other, the entropy must stay constant because the evolution of the quantum state 
is unitary. As the relation between entropy and v is fixed for ME-packets, the ME- 
packet form is not preserved by the evolution (the entropy ceases to be maximal). 
This is similar for Gaussian-packet form or for coherent-state form. 

For general potentials, there will be two types of corrections to the dynamics of 
the averages: terms containing the variances and terms containing h. To see these 
corrections, let us calculate time derivatives for the Hamiltonian ( TlTl) with potential 
( 1291) . The Heisenberg-picture equations of motion give again 

dq 1 

at m 

so that Eq. (1301) is valid. The other equation, 

can be applied iteratively as in the classical case so that all time derivatives of p can 
be obtained. Thus, 

d i = -V 1 -V 2q -^ g ^^ + r 5 , (74) 
at 2 b 

and 

d 2 p V 2 V 3 V 4 2 2 

— = p - —{qp + pq) - —{q p + qpq + pq ) + r 5 . 

at z m 2m bm 

This differs from the classical equation only by factor ordering. We can use the 

commutator [q,p] = ih to simplify the last term, 

72 



d 2 p V 2 V 3 . . V 4 



lr > -—v - i^{qp + pq) - 7r-qpq + r s • (75) 

at z m 2m 2m 



Similarly, 



d 3 p V 3 2 V 4 V X V 2 V X V 3 + Vj 3V 2 V 3 + 2 

-T7g = 2 P zpqp + + q + ~ q 

at 6 m A m z m m 2m 

AV 2 V 4 + 3Vi 3 5V 3 V, 4 Vl 5 

h t. -q + — — q + 77— q + r 5 , (76) 



bm ' 12m 12m 



and 



-771 = 3 P + P + JT~2 ( qP + Pq ) + ^ ^ 

ctt 4 m 6 m z 2m z 2m z 

3V 3 Va , q o, 3V? o o / \ 

+ "2^ (g P + ^ + 4m? + rs • (7?) 
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Next, we calculate quantum averages with the help of formula (ITT]) . The quantum 
averages of the monomials that are linear in one of variables q or p can differ from 
their classical counterparts only by terms that are of the first order in \jv and purely 
imaginary. For example, 

(qp) = QP + j , 



or 



(q 3 p) = Q 3 P + 3QPAQ 2 + 3i 



Q 2 AQAP AQ 3 AP 

h 61 



V V 

These corrections clearly cancel for all symmetric factor orderings. The first term 
in which a second-order correction occurs is q 2 p 2 and we obtain for it: 

/ 2 v / 2 2v AQ 2 AP 2 

(pq p) = (q p ) class + 2 . 

v A 

The equations ( 174l) -( l77l) do not contain any such terms and so their averages 
coincide exactly with the classical equations (I35l) - (l38l) . The terms q 2 p 2 with different 
factor orderings occur in the fifth time derivative of p and have the form 



3^4 



2m 2 



J) 

q 3 p + pq 3 ,— 
2m 



+ 



V3V4 
2m? 



1 o o 



ih 



^(21pq 2 p-llh 2 ) 



The average of the resulting term in the fifth time derivative of p is 



V3V4 
2m? 



2\Q 2 P 2 + 21P 2 AQ 2 + 2ig 2 AP 2 + 21AQ 2 AP 2 - 



h 1 



If we express h as 2AQAP/ u, we can write the last two terms in the parentheses as 

AQ 2 AP 2 (21 - 4 

A similar term appears in the third time derivative of p, if we allow V5 7^ in the 
expansion ( l29l) : 



{q p + pq ),— 
2m 



12m 

which contributes to d 3 P/dt 3 by 

2m 2 



' (2pq 2 p 



Am 



h 2 ) 



(q V ) class + 



4AQ 2 AP 2 



Again, the correction is of the second order in u^ 1 . 

We can conclude. The quantum equations begin to differ from the classical one's 
only in the higher order terms in V or in the higher time derivatives and the correc- 
tion is of the second order in \ jv. This seems to be very satisfactory: our quantum 
model reproduces the classical dynamic very well. Moreover, Eq. ( 1691) shows that 
Gaussian wave packets are special cases of ME packets with v = 1. Thus, they 
approximate classical trajectories less accurately than ME packets with large v. 
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4 Classical limit 



At some places of the paper, it is written that v ^> 1 is the classical regime. Let us 
now look to see if our equations give some support to this statement. 

Let us consider averages of powers of q and p. If we expand such an average in 
powers of Planck constant then the leading term can be calculated with help of the 
formula (flOl) from quantum partition function ( J60l) and from relations ( JoTl) and ( l62l 
between the Lagrange multipliers and the averages and variances of q and p. This 
has been explained after the proof of Lemma 2. 

The quantum partition function (1601) differs from its classical counterpart ([6]) by 
the denominator sinh( hy/X^Xl). If 



h^X, < 1 , (78) 

we can write 

The leading term in the partition function then is 

71 1 ( X\ \l 



~ h v^A^ 6XP V4A^ 4A 4 

where h = 2irh. Comparing this with formula (JSD shows that the two expressions 
are identical, if we set 

v = h . 

We can say that quantum mechanics gives us the value of v. Next, we have to 
express condition (ITHj) in terms of the averages and variances. Equations ([52]) imply 

fr\A3A 4 = r In V -^~ ■ 
2 v — 1 

Hence, condition ( 1781 is equivalent to 

v > 1 . (79) 

The expression 

f + 1 

— In 

2 v-\ 

that appears on the right-hand sides of Eqs. AMI) and (162!) satisfies 

lim — In — -t— ■ = 1 . 

i^=oo 2 z/ — 1 

Hence, the leading terms in these equations coincide with Eqs. ([Tj) and jSJ). 
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Our result can be formulated as follows. The time evolution of classical and 
quantum ME packets with the same initial values of averages and variances defines 
the averages as time functions. These time functions coincide for the two theories 
in the limit v = oo. Hence, in our approach, this is the classical limit. It is very 
different from the usual assumption that the classical limit must yield the variances 
as small as possible. One also often requires that commutators of observables van- 
ish in classical limit. This is however only motivated by the assumption that all 
basic quantum properties are single values of observables. Within our interpreta- 
tion, this assumption is rejected and if classical observables are related to quantum 
operators then only by being average values of the operators in prepared states. All 
such averages are defined by the preparation and do exist simultaneously, indepen- 
dently of whether the operators commute or not. For example, Q and P are such 
simultaneously existing variables for ME packets. 

Let us compare the present paper notion of classical limit with a modern textbook 
version such as Ch. 14 of [22] . Both approaches define the classical limit of a quantum 
state as a classical ensemble described by a fuzzy distribution function and calculate 
time evolutions of averages in the states. However, in the textbook, any quantum 
system, even not macroscopic, and any state, even pure, are allowed (pure states 
preferred as they have smaller uncertainties). Hence, our notion is much narrower: 
we consider only macroscopic quantum systems and only some of their maximum 
entropy states. This has obvious physical reasons explained in the first two sections. 

5 Conclusion 

The paper describes a quite general construction of quantum states that model 
important properties of classical-mechanical states. To achieve that, one often- 
assumed classical property has to be abandoned: the completely sharp trajectory of 
all mechanical systems. The sharp trajectory is considered here only as an ideal limit 
allowed by classical mechanics. There is however nothing in nature that corresponds 
to it. This is in agreement both with practical observations and with theoretical 
idea that the correct underlying theory is quantum mechanics. Hence, the way 
for statistical methods highlighted in Ref. pE] is free. The key concept turned out 
to be physical conditions equivalent to preparation process in quantum mechanics. 
The paper transfers it into classical mechanics, where it generalises the old notion 
of initial data. Entropy is defined by the physical conditions independently of the 
state of any observer's mind. 

Classical mechanics allows not only sharp, but also fuzzy trajectories and the 
comparison of some classical and quantum fuzzy trajectories shows a very good 
match. The fuzzy states chosen here are the so-called ME packets. Their fuzziness 
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is described by the quantity v = 2AQAP/ h~. The entropy of an ME packet depends 
only on v and is an increasing function of it. The larger v is, the better the quantum 
and the classical evolutions of average values have been shown to agree. Thus, the 
classical regime is neither AQ = AP = (absolutely sharp trajectory) nor v — 1 
(minimum quantum uncertainty). This is the most important result of the paper. 

Unlike internal classical properties the external ones such as coordinates and mo- 
menta are well manipulable so that conditions exist allowing ME packets to have all 
variances AQ and AP from a broad range. There is no reason why the variances 
had always to be small other than the incorrect assumption that all real mechani- 
cal trajectories are absolutely sharp. Hence, our account of classical properties as 
statistical properties of macroscopic quantum systems that started in [4] can be 
considered as concluded. 
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